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In hybrid inflation, initial field values leading to sufficiently long inflation were thought be 
fine-tuned in a narrow band along the inflationary valley. A re-analysis of this problem has 
shown that there exists a non negligible proportion of successful initial conditions exterior to 
the valley, organized in a complex structure with fractal boundaries, and whose origin has been 
£N| ' explained. Their existence in a large part of the parameter space has been demonstrated using 

a bayesian Monte-Carlo-Markov-Chain (MCMC) method, and natural bounds on potential 
parameters have been established. Moreover, these results are shown to be valid not only for 
the original hybrid model, but also for other hybrid realizations in various frameworks. 
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1 Introduction and original hybrid model 

O , 

Q" 1 Recent observations of the Cosmic Microwave Background (CMB) and its anisotropics have 

provided strong arguments in favour of a phase of accelerated expansion in the early universe. 
If the simpler way to realize this inflation era is to assume the universe initially filled with 
an unique homogeneous scalar field slowly rolling along its potential, many other realizations 
have been proposed §. In hybrid models, an inflaton field is coupled to an auxiliary field, and 
accelerated expansion usually occurs when inflaton slow-rolls along a nearly flat valley of the 
potential and ends abruptly due to a tachyonic instability. The advantages of hybrid models are 
that the energy scale of inflation can be low and initial field values do not need to be larger than 
the Planck mass. Moreover, the model can be embedded in some high energy frameworks like 
supersymmet rv, supergravity and grand unified theories. A generic prediction of the original 
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hybrid modeV^ 1 is a blue spectrum of scalar perturbations, and thus it is strongly disfavoured 
by WMAP constraints on the scalar spectral index^. It remains nevertheless a good toy model 
for other hybrid realizations in various frameworks, whos e dynamics is similar. 
^ \ The potential for the original hybrid model reads^El 
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in which <j) is the inflaton, tp is the auxiliary Higgs-type field, and M, fi, 4> c are three mass 
parameters. Inflation occurs in the false-vacuum along the nearly flat valley tp = 0. A tachyonic 
instability appears when inflaton reaches <j> = cp c . From this point, the trajectory falls through 
one of the global minima of the potential d> = 0, ib = ±M whereas tachyonic preheating occurs 

It is a natural question to ask how the fields have to be tuned initially along the inflationary 
valley in order to generate more than around 60 e- folds of inflation inside the valley. If trajectories 
are initially displaced slightly in the transverse direction, they perform damped oscillations along 
the valley before the slow-roll regime engages the realization of a large number of e-folds. On the 
contrary, if the initial value of the auxiliary field is too large, the damping can not proceed 
sufficiently quicklyand slow-roll regime do not begin before the critical point of instability is 
reached. Tetradis^ and afterward Mendes and Liddled determined that the successful region 
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in the initial field space was a very narrow band along the valley such that the initial value 
of the auxiliary field had to be fine-tuned compared to initial inflaton values. Howev er, som e 
inconsistencies can be pointed out and have lead to the reanalysis of the problem inEMIOl _ 
Both El and^ observed some unexplained successful initial conditions outside the valley, who 
seem organised in some structures in^, but apparently isolated and very subdominant in the 
higher-resolution analysis of^ . These previous studies were also restricted to some similar sets 
of potential parameters. A full quantitative study of the problem of initial field values in the 
whole parameter space was still lacking. 

By integrating numerically the classical 2-field dynamics, the behaviour of such successful 
trajectories starting outside the valley can be understood. The relative area that these points 
can occupy can be evaluated for some sets of parameters, and their dominance in the whole space 
of potential parameters, initial field-values and initial velocities can be probed using a statistical 
Monte-Carlo-Markov-Chain method. Our results indicate that the model do not suffer of fine- 
tuning any more because successful trajectories start more probably outside the valley. It is 
also shown that initial velocities do not affect the probability to generate sufficient inflation. 
Finally, some bounds on potential parameters can be established with the only requirement of 
sufficiently long period of inflation and the absence of some fine-tuning of initial conditions. 
The important questions of the robustness of these observations for other hybrid realizations, 
as well as the effects of quantum fluctuations on such trajectories can be pointed out and will 
be addressed briefly in the conclusion. 



2 Space of initial conditions 



2. 1 Grids of initial field values for fixed parameters 

High-resolution grids of initial field values have been plotted by integrating numerically the clas- 
sical 2-field dynamics. In a flat FLRW universe, dynamics is governed by Friedmann-Lemaitre 
equations 
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as well as Klein-Gordon equations in an expanding universe 
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where a dot denotes derivative with respect to cosmic time, m p \ is the Planck mass iff = a/a 
and a is the scale factor. In a first step, initial velocities are assumed to be vanishing. White 
regions in the grid represented in figure [1] correspond to initial conditions leading to more than 
60 e-folds of inflation. As expected, a band of fine-tuned initial field values is found along the 
valley, as well as successful points outside the valley. But instead of being isolated, as it was 
suggested in^ , they form a complex structure of thin lines and crescents, and with successive 
zooming it is observed that the structure form a connected set. Moreover, the fractal properties 
of these successful initial conditions have been studied El. The box-counting dimension of the 
boundaries has been determined numerically and found to be non-integer. Nevertheless, the 
structure itself exhibits a non-fractal finite area. These properties are similar to the well known 
Mandelbrot set whose area is finite with infinite fractal boundaries. 

Finally, the relative area covered by successful points have been quantified. These points are 
found to cover up to 20% of the space of initial field values, depending on potential parameter 
sets 1^1. Therefore, the amount of fine-tuning necessary to have successful inflation is strongly 
reduced for these particular sets of parameters. 



b M p i is used for the reduced Planck mass m p i/\/87r 
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Figure 1: grid (2040x2040 points) of initial field values producing more/less than 60 e- folds of inflation (white/black 
regions). The potential parameters are set to M = 0.03 m p i, = 636 m p i, <f> c — \/2 x 10~ 2 m p i. 

2.2 Successful trajectories initially outside the valley 

The typical behaviour of successful trajectories starting outside the inflationary valley (see figE]) 
has been studied more precisely. Three successive phases have been identified. First, a fast- 
rolling phase pushing the trajectories through the bottom of the potential, without generating 
a significant number of e-folds. Then, after some rebonds on the sides of the potential, these 
trajectories become oriented along the valley, and climb it. Arriving at an extremum point with 
a quasi- vanishing velocity |, they start to slow roll back along the valley while producing a large 
number of e-folds. 

2.3 Statistical MCMC analysis 

Grids of initial conditions demonstrates that a severe fine-tuning of initial field values is not 
required in hybrid inflation for particular sets of potential parameters. However, the full initial 
fields and potential parameter space including initial velocities is 7-th dimensional @ and thus is 
hardly explored with such grids. For this reason, a statistical bayesian MCMC analysis of this 
space has been performed. Priors are chosen to be flat for initial field values, initial velocities, 
and for the logarithm of the potential parameters, in order to avoid fine-tuning of the initial 
fields and to not favour any precise scale for the potential parameters. Only the sub-Planckian 
regime has been considered, but it was shown in^ that for super-Planckian trajectories, the 
realization of at least 60 e-folds of inflation becomes generic. More details about the MCMC 
method can be found in^. Only the main results will be discussed below. 

c This extremum point is located inside the fine-tuned successful band, thus there exists a correspondence 
between each initial value inside and exterior to the valley. 

d Let notice that A has not been considered since it only normalizes the potential without influencing the 
dynamics. 
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Figure 2: Typical behaviour of successful trajectory starting outside the inflationary valley, first falling towards 
the bottom of the potential, then oscillating and climbing the inflationary valley (tp — direction) instead of 
being trapped by one of the global minima of the potential (represented by crosses). Inflation occurs when it 

slow-rolls back along the valley. 



The first result concerns the posterior probability density distribution of tf)\ marginalized 
over the rest of the space (see figure [3]) . It exhibits a maximum at ip\ = corresponding to 
successful trajectories starting inside the valley. However, it is important to remark that the 
probability distribution is widely spread over large values of ipi, so that the probability to have 
a successful trajectory starting exterior to the valley is dominant over the probability to start 
along it. It can therefore be concluded that in the original hybrid model, there is no need of 
fine-tuned initial conditions along the valley to generate sufficient inflation. Concerning initial 
velocities, posterior probability distributions are flat, and they do not play a role in the game 
of finding trajectories leading to many e-folds of inflation. Finally, from the only requirement of 
having more than 60 e-folds of inflation, natural bounds on the potential parameters (ft c and \x 
can be established, 

<p c < 0.004 m p i 95%C.L., fi > 1.7 m p i 95%C.L.. (4) 

From the first bound, inflaton value at instability point is shown to be smaller than the reduced 
planck mass, and the last 60 e-folds are sure to be generated in the sub-planckian regime inside 
the valley. The second bound is linked to the non-existence of a small field phase of inflation 
along the valley for small values of /x. This regime was put in evidence in^. 
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Figure 3: Posterior probability density distribution for ipi (in m p i units) with the requirement of at least 60 e-folds 
generated, marginalized over initial values of <f>, potential parameters and initial field velocities. 



3 Conclusion and discussion 



By integrating the classical 2-field dynamics, it has been shown that the original hybrid model 

does not require fine-tuned initial field values in order to generate a sufficient number of e-folds 

to solve the standard cosmological problems. Actually, trajectories do not need to be initially 

tuned inside a narrow band along the inflationary valley. They can start exterior to the valley, 

and after a fast-roll phase and some rebonds on the sides of the potential, they join it and enter 

in the slow-roll regime. A bayesian MCMC method has been used to show that this observation 

is valid, not only for some specific potential parameter sets, but in a large part of the 7-D space of 

initial field values, initial velocities and potential parameters. Initial velocities do not influence 

the probability to generate sufficient inflation, and bounds on potential parameters have been 

deduced from the only requirement of a sufficiently long period of inflation. 

If the analysis of the successful initial conditions has been conducted essentially for ori gina l 
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hybrid model, the results remain generic for other hybrid realizations in various framework 1 ^ 121 . 
A large part of the space of initial fie ld valu es (up to 80%) has been found to be successful for 
super symmetric smooth^ and shifted^2H31 inflation, both in their SUSY and SUGRA versions, 
as well as for radion assisted gauge inflation^]. The MCMC analysis have been conducted for 
F-term hybrid inflation^^ in supergravity, with similar results. In particular, a bound on the 
only parameter M of the model has been established, 



log(M/M p i) < -1.33 95% C.L., (5) 

comparable to the upper bound obtained from cosmic strings formation in^^l 

One could argue that there is a remaining problem in the sense that successful initial con- 
ditions are subdominant (they do not exceed 25% of the space of initial field values for original 
hybrid model) compared to unsuccessful trajectories. However, it is without noticing that suc- 
cessful points outside the valley form a complex fractal structure covering roughly all parts of 
the space of initial conditions. In the context of self-reproducing universe, due to quantum fluc- 
tuations along trajectories, each initial patch of the universe could contain at least one region 
for which inflation will undergo while in other regions inflation will not occur. At the end, the 
volume occupied by the regions in which inflation occurred is largely dominant. This effect 
should be studied more in details in a future work. 
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